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Abstract 

- Context The frequency, size, and insertion angle of 
primary branches are important determinants of wood 
quality and can be significantly influenced by silvicultural 
activities. 

- Aims This study quantified the long-term influence of 
early re-spacing on the branching characteristics of mature 
Sitka spruce [Picea sitchensis (Bong.) Carr.]| trees growing 
in Northern Ireland, UK. The primary aim was to 
investigate whether any residual effect of stand density on 
branch attributes remained once the effect of changes in tree 
size variables had been considered, while a secondary 
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objective was to test the performance of existing Sitka 
spruce branch models using the current dataset. 

- Methods Re-spacing treatments had corresponding stand 
densities of 2,858, 1,452, 725, 477, and 320 stems ha |. 
Twenty-four trees were sampled when the stand was 
57 years old and branch frequency, size, and insertion 
angle were recorded for model development. 

* Results Maximum branch diameter, insertion angle, and 
branch frequency were significantly influenced by re-spacing, 
while no effect was found for relative branch diameter 
distribution. Residual re-spacing effects were most noticeable 
on branch size, with only small differences between treat- 
ments for branch frequency and insertion angle. Existing 
models performed well despite the wider range of stand 
densities examined in the present study. 

- Conclusion The results indicate that early re-spacing from 
1.9 m* to wider than 2.6 m* will result in branch attributes 
that are detrimental to Sitka spruce sawn timber quality. 


Keywords Picea sitchensis - Wood quality - Early 
re-spacing - Maximum branch size - Number of branches - 
Branch angle - Relative branch size 


1 Introduction 


Sitka spruce [Picea sitchensis (Bong.) Carr.] is the most 
commonly planted conifer species in the UK (Forestry 
Commission 2010) and Ireland (Joyce and OCarroll 2002) 
and is the main commercial species processed by the 
sawmilling and panel board industries in both these 
countries. Due to a significant increase in planting during 
the period spanning the 1950s to the 1970s, timber 
production is expected to gradually increase over the next 
15-20 years in line with the current UK forecast (Halsall 
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et al. 2006). Since Sitka spruce will account for over 80% 
of this increase, it will be necessary to increase the 
proportion of locally produced timber used in the larger 
and potentially higher-value construction market, as exist- 
ing markets for nonstructural uses such as pallets, fencing, 
and packaging have become saturated. For this to occur, a 
greater proportion of Sitka spruce timber will have to meet 
the quality requirements demanded by end-users, which 
include strength class, as defined in Europe by EN338 
(CEN 2003), dimensional stability, and visual appearance. 

The quality of sawn timber is determined by the 
properties of the logs from which it is produced. For 
structural end uses, important wood properties include 
wood density, microfibril angle, and the size and frequency 
of knots, as these affect wood stiffness (MOE) and strength 
(MOR; Forest Products Laboratory 2010) and hence, the 
structural grade outturn under both visual and machine 
strength grading systems (Maun 1992; Zhang et al. 2002). 
In addition, problems can also occur when using exces- 
sively knotty timber in timber frame construction due to 
nail fouling. The characteristics of knots are, in turn, 
determined by branch attributes such as frequency, size, 
insertion angle, and status. Since these attributes are 
strongly influenced by growth characteristics and competi- 
tion, there is potential for these to be modified by forest 
managers in order to improve the yield and performance 
of structural timber (Brazier 1977; Macdonald and 
Hubert 2002). 

However, there are also concerns about the effects on 
timber quality of changes in UK silvicultural practices that 
were introduced in the 1960s and continued through to the 
1980s, such as a move to wider initial spacings (Wardle 
1967), partly to reduce establishment costs but also to 
improve the wind stability of stands, and increased planting 
of upland areas (Worrell 1987). These silvicultural trends 
may have had a detrimental effect on future Sitka spruce 
timber quality because increased tree growth is associated 
with larger knots, reduced wood density, and poorer stem 
form (Brazier 1970, 1977). For example, a recent study 
investigating the effects of an early re-spacing treatment 
on the physical and mechanical properties of Sitka spruce 
structural timber found that wide re-spacing resulted in 
a significant reduction in both MOE and MOR (Moore 
et al. 2009). 

The wide range of silvicultural treatments applied by 
forest managers due to the competing objectives of multi- 
use forestry can influence timber quality both positively 
and negatively (Macdonald et al. 2010). To account for this 
variability, a modelling approach has evolved that links 
wood quality with tree growth in order to quantify 
predictions of end-use properties under different silvicul- 
tural scenarios (Houllier et al. 1995). Specifically, this 
approach enables wood quality attributes such as density 
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(e.g. Guilley et al. 2004; Schneider et al. 2008; Gardiner 
et al. 2011), stem taper, (e.g. Tasissa and Burkhart 1998; 
Garber and Maguire 2003; Fonweban et al. 2011), and 
branch attributes (e.g. Maguire et al. 1991; Colin and 
Houllier 1992; Weiskittel et al. 2010) to be modelled as 
functions of tree diameter and height growth as influenced 
by silvicultural practices. 

Using such a modelling approach, Achim et al. (2006) 
studied the changes in predicted branch characteristics of 
Sitka spruce following thinning through the effect on tree 
size variables (i.e. diameter and height). The authors found 
that while thinning had no influence on branch angle, it 
significantly increased branch frequency and size. However, 
this analysis was limited to the range of initial stand densities 
commonly used in the UK (e.g. approximately 1,750-— 
3,500 stems ha‘; Achim et al. 2006), therefore it is 
uncertain whether the models will give accurate predictions 
for trees grown outside the range of conditions for which 
they were calibrated, even though such extrapolations are 
commonly made (Macdonald et al. 2010). In addition, 
Achim et al. (2006) did not explicitly consider the effect of 
stand density over and above the effect on tree growth (i.e. 
the marginal or additional effect of thinning treatment after 
the influence of tree size variables had been accounted for). 

Therefore, the overall goal of this study was to quantify 
the influence of a wide range of re-spacing intensities on 
Sitka spruce branch attributes, in order to test the 
assumption that they can be explicitly modelled by linking 
them to tree growth. More specifically, we used a long- 
established early re-spacing trial with a wide range of stand 
densities in order to test the hypothesis that there is no 
influence of re-spacing intensity on branch angle, frequency, 
and maximum as well as relative branch size in Sitka spruce 
over and above the effect that can be modelled using changes 
in tree size variables such as diameter at breast height or 
relative crown dimensions. A secondary objective was to 
establish whether existing models (1.e. Achim et al. 2006) 
can be successfully extrapolated beyond the range of stand 
densities for which they were originally developed. 


2 Methods 
2.1 Experimental data 


Twenty-four trees were sampled from a re-spacing exper- 
iment at Baronscourt, in Co. Tyrone, Northern Ireland 
(latitude 54° 41'N, longitude 7°26'W, 140 m elevation). 
The 4-ha experimental site, with a yield class of 
20 m° ha ' year ', was planted in 1949 at an initial spacing 
of 1.83<1.83 m, and re-spaced in 1960. Five re-spacing 
treatments, which included the initial spacing, were applied 
to 25 plots of approximately 0.15 ha in size arranged in a 
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5x5 Latin square design. In addition to the initial spacing, 
re-spacing treatments were nominally 1.83x3.66, 3.66 
3.66, 3.66<5.49, and 5.49x5.49 m, and the corresponding 
measured stand densities were 2,858, 1,452, 725, 477, and 
320 stems ha -, respectively. 

For brevity, we refer to the treatments by the equivalent 
Square spacing used by Kilpatrick et al. (1981), Le. 
approximately 1.9, 2.6, 3.7, 4.6, and 5.6 m’*. The four 
wider re-spacing levels were achieved by systematically 
removing alternate rows and/or trees within a row. The 
experiment is described in more detail in Jack (1971) and 
growth and yield data for the first 31 years are presented in 
Kilpatrick et al. (1981). Due to significant mortality, the 
stand densities in the 1.9 and 2.6 m* treatments at the time 
of sampling were significantly lower than in 1960 (Table 1). 
Due to past wind damage, one of the five plots in the 
4.6 m* spacing treatment was omitted as the growing space 
available to the trees was no longer representative of this 
treatment. 

In 2006, when the stand was 57 years old, the diameter 
at breast height (DBH) of every live tree in all plots was 
recorded. One dominant or co-dominant tree from each plot 
was randomly selected from the upper 50th percentile of the 
DBH distribution for destructive sampling. The DBH of 
these trees ranged from 33 to 76 cm, and total height (HT, 
measured once the tree had been felled) ranged from 29.2 to 
38.6 m (Table 1). Detailed branch attributes were measured 
using the approach described in Colin and Houllier (1991) 
and Achim et al. (2006). 

In Sitka spruce, branches tend to form either in a distinct 
cluster in the upper portion of each annual shoot—more 
correctly termed ‘pseudo-whorls’ in the botanical literature 
(e.g. Heuret et al. 2006)—or along the length of each 
annual shoot. Here, we refer to these as nodal and 
internodal branches, respectively (Barthélémy and Caraglio 
2007). On each annual shoot (i.e. the continuous elongation 


Table 1 Stand-level characteris- 


tics in 2006 for each re-spacing Attribute 
treatment and characteristics of the 
24 sample trees on which branch 
attributes were measured 
Stand 


Stand density (stems ha ') 
Mortality (%) 


Where iate, standard 
re appropriate, standar DBH (cm) 


errors of the mean are given in 
parentheses. Standard errors for Height (m) 
stand-level characteristics are Sample tree 
based on data from plots within N 

a re-spacing treatment 


DBH diameter at breast height ee 
(cm), HT total tree height (m), Height (m) 
LLPW height (m) of the lowest LLPW (m) 
live pseudo-whorl of branches, LLB (m) 


LLB \owest live branch height (m) 


oD 


flush of the stem during a single growing season), the 
distance (DIST, metre) from the shoot base to the stem apex 
was measured and the annual shoot length (i.e. the 
internodal length, NLEN) was by calculated by subtracting 
DIST values for two successive annual shoots. On every 
fourth annual shoot, the position (i.e. nodal or internodal), 
status (i.e. alive or dead), vertical insertion angle (i.e. the 
angle of insertion relative to the upright stem, VBA), and 
horizontal and vertical diameters (near branch base, avoid- 
ing basal swelling) of each branch greater than 5 mm in 
diameter were recorded. Branch diameter (BD) was 
calculated as the arithmetic mean of the horizontal and 
vertical diameters. The relative branch diameter (BDppfy) of 
each branch was determined as the diameter relative to the 
maximum branch diameter in each annual shoot (BD/ 
BDwmax). Crown length (CL) was defined as the distance 
(m) between the stem apex and the insertion point of the 
lowest living branch (i.e. the lowest branch containing live 
foliar tissue), and crown ratio (CR) was calculated as CL/ 
HT. Overall, a total of 2,528 branches were assessed on 295 
annual shoots (12.3+0.7 annual shoots per tree). 


2.2 Data analysis 


Because the data had a hierarchical structure with branches 
clustered within annual shoots, that were themselves 
located within trees, a mixed-modelling approach similar 
to that employed by Hein et al. (2008a) was adopted to 
ensure that appropriate estimates of parameter standard 
errors were obtained and tests of parameter significance 
were valid (Pinheiro and Bates 2009). A linear mixed- 
effects modelling approach was used to analyse data on 
maximum branch diameter, relative branch diameter, and 
branch insertion angle. Because the number of branches per 
annual shoot represents count data and the mean counts per 
branch type were generally low (1.e. <5), a generalised 


Re-spacing treatment (m7) 


2.6 a7 4.6 5.7 

1,134 (71) 922 (150) —-581 (45) 435 (12) 310 (19) 
60 37 20 9 3 

29.1 33.5 41.3 46.5 512 

32.8 34.5 34.9 36.2 35.8 

5 5 4 5 5 

38.9 (5.7) 36.9 (5.1) 403(6.9)  55.6(5.7) 65.9 (9.1) 
32.5 (2.6) 34.7(1.2)  34.6(2.1) 32.5(2.6) 37.0 (1.5) 
23.1 (1.8) 2436) 24506) B38) D105) 
57) 23.6 (1.8)  23.4(2.0) 223(1.8) 20.3 (1.1) 
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linear mixed-effects modelling approach, which assumes 
that counts follow a Poisson distribution (McCullagh and 
Nelder 1989), was used to model the number of branches. 
In the models for branch number and maximum branch 
diameter, random effects were included for tree, while 
random effects for tree and annual shoot were included in 
the models for relative branch diameter and branch 
insertion angle. 

For each branch characteristic of interest, the final choice 
of explanatory variables to be included in the fixed part of 
the models were determined by comparing the results from 
various initial model fits using Akaike’s information 
criterion (Akaike 1974) with a difference of 10 units being 
considered significant (Burnham and Anderson 2002, p. 
71). In addition, data on the branch and annual shoot 
characteristics were treated as repeated measures along the 
stem or along a gradient within the annual shoot. 
Preliminary analysis of the data indicated there was a 
significant longitudinal correlation, and a continuous 
autoregressive model was included to account for this 
(Pinheiro and Bates 2009). The significance of both the 
random effects terms and the correlation function in each 
model were determined using likelihood ratio tests. Only 
those parameters that were significant (p<0.05) were 
retained in the final models, except where interaction terms 
were included in the models, in which case all fixed-effects 
parameters are reported. Multiple comparisons between 
re-spacing treatments were carried out using Tukey’s 
adjustment. All data analyses were performed using the 
functions contained in the Ime4, nlme, and multcomp 
libraries of the R_ statistical software programming 
environment (R Development Core Team 2011), except 
the maximum branch diameter model which was fitted in 
Statistical Analysis System (SAS) using the PROC MIXED 
procedure (Littel et al. 2002). In the latter case, SAS was 
used rather than R as the latter currently lacks the facility to 
carry out multiple comparisons when significant model 
interaction terms that include the factor of interest are 
encountered. 

Initially, all models were developed without tree or 
crown size variables in the fixed-effects predictors, as these 
are affected by stand density. First, base models were 
developed for each branch characteristic of interest, and 
included re-spacing as a categorical factor in the fixed 
effects. The 1.9 m? re-spacing treatment was used as the 
reference treatment in all models. Next, tree or crown size 
variables were added to the fixed effects of each model in 
order to determine the marginal effect of re-spacing 
intensity on branch characteristics once any effects attrib- 
utable to changes in tree size and crown dimensions had 
been accounted for. Several tree- and crown-level variables 
were screened, including DBH, CL, CR, and HT/DBH 
ratio. Although there were minor differences in the number 
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of trees in each plot, preliminary analysis (results not 
shown) indicated that the inclusion of re-spacing treatment 
as a categorical factor was a more effective approach than 
using a continuous measure of stand density. Because data 
from both live and dead branches along the stem were used 
in developing the final equations, any effect of branch 
status on branch attributes was evaluated by adding an 
indicator variable (1 if live branch, 0 if dead branch) to 
each model. 

For comparative purposes, the model forms and param- 
eter estimates developed by Achim et al. (2006) were also 
applied to this dataset in order to predict branch number 
and insertion angle. Since the model for branch diameter 
developed by Achim et al. (2006) predicts the mean branch 
diameter in an annual shoot, it could not be compared to the 
model of maximum branch diameter developed in this 
study. 


2.2.1 Branch number 


Separate models were developed to predict the number of 
nodal (NBy) and internodal branches (NB,) in an annual 
shoot, as well as the total number of branches (NB7). 
Following the approach taken by Hein et al. (2008a), 
branch number in each annual shoot for each branch type 
was related to both the annual shoot length and the annual 
shoot height in the stem: 


NBy = ap + a3NLEN+a.NHT-+a3DBH-+a4SP) 6+a5SP3 7 
+a6SP4.6+a7SPs5.7+0tree 


(1a) 
NB, = bp + 6, NLEN+6.NHT + 63DBH + 64SP> 6 


+ bsSP3.7 + beSP4.6 + 57SP5.7 + Otree (1b) 


NBr = co + cj} NLEN + @NHT + c3DBH + c4SP 6 
+ csSP3.7 + c6SP4.6 + C7SP5.7+0tree (1c) 


where NLEN is the annual shoot length (m), NHT is the 
height in the stem (m) of each annual shoot (i.e. HT — DIST), 
DBH is tree diameter at breast height (cm), a ...a7, b,...b7 
and co...c7 are parameters to be estimated from the data and 
Otree 18 a random effect for each tree. SP> 6, SP37, SP46, and 
SP; ¢ are indicator variables for the 2.6, 3.7, 4.6, and 5.6 m* 
re-spacing treatments, respectively. 

The number of nodal and internodal branches in an 
annual shoot was also predicted for the model-fitting 
dataset using the following equation developed by Achim 
et al. (2006): 


NBy s = iy, NLEN*¥! e(-?vs PST) (2) 
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where NLEN and DIST are as previously defined and iy ;, 
ay and by were modelled as functions of the tree-level 
variables HT and DBH for both nodal and internodal 
branches, 1.e.: 


in = 5.71 —0.0084HT 77; = 12.49 — 0.18 


HT 


HT 
by = —0.089 + 0.072DBH + 0.0055 —— __ by; 


Spur 01 = 0.23 


where DBH* is in metres to ensure correspondence of units. 
2.2.2 Maximum branch diameter 


The diameter of the thickest branch in each annual shoot 
was modelled using the following linear equation: 


BDmax = dp + d,DIST + dy In(DIST) + d3CR 
+:049P36 1 G59P3.9 + dooPae + G79F 56 
+ de(SP26 x DIST) + do(SP3.7 x DIST) 
+ dio(SP46 X DIST) 
+ d\,(SPs5.6 X DIST)+6 tec (3) 


where BDyax iS the maximum branch diameter (mm), 
DIST is the distance of each annual shoot from the stem 
apex (m), In(DIST) denotes the natural logarithm of DIST, 
SPxDIST is an interaction term, do...d;, are parameters to 
be estimated from the data and dt. represents the tree-level 
random effect. 


2.2.3 Relative branch diameter 


The relative diameter of all branches in an annual shoot is 
an important attribute because it describes the size variation 
of branches on a stem and allows the diameter of any 
branch in an annual shoot to be estimated if the number of 
branches and the maximum branch diameter are known. 
The relative diameter of all branches in an annual shoot, 
except the largest, was modelled using a modified version 
of the linear equation presented in Hein et al. (2008a): 


BDrerz =p e,RANK + €@9 In(RANK) + esNBr 
+ e4NHT + esN + e6SP2.6 + e7SP37 
5 €g SP 4.6 ie €oSP5 6 ae Ogee = Onbot (4) 


where BDrey is branch relative diameter within an annual 
shoot (BD/BDyax), RANK is the branch rank within an 
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annual shoot (e.g. 1 for BDyax, 2 for the next largest branch, 
etc.), In(RANK) the natural logarithm of RANK, NB, is the 
total number of branches within the annual shoot, N is an 
indicator for branch type (1 if nodal branch, 0 if internodal 
branch), and eg...@¢9 are parameters to be estimated from the 
data. Oghoot 18 the random effect for each annual shoot. 


2.2.4 Branch insertion angle 


The angle of insertion of a branch relative to the vertical 
stem was modelled using the following linear equation: 


VBA = go + gi BD + g2DIST + g3 In(DIST) 
+ g4(BD x DIST) + gsN + gol + g7SP2.6 
+ geSP3.7 + goSP4.6 + 2105P5.6 
+ Otree FO shoot (5) 


where VBA is the vertical branch angle (°), BD DIST is an 
interaction term, Z is an indicator variable for branch status 
(1 if live branch, 0 if dead branch), go...g19 are parameters 
to be estimated from the data, and all other variables are as 
defined previously. 

VBA of nodal (VBAj) and internodal (VBA,) branches 
was also predicted for the model-fitting dataset using the 
following equations developed by Achim et al. (2006): 


__ 0.0028 


VBAy = 72.87e7 1-2 (6a) 





VBA, = 85.6e7 039-2 (6b) 
where Z is the relative distance of each annual shoot from 
the stem apex (DIST/HT). 


2.2.5 Evaluation of model performance 


The performance of the fixed part of each of the models 
was examined by calculating the mean error (£), mean 
absolute error (|E]), and mean percentage error (E%) using 
the model-fitting dataset. E and |E| were calculated using 
the equations presented in Hein et al. (2008a), while E% 
was calculated using the following equation: 


100 ys |vi 9, | 


where y; is the observed value, ); the predicted value, and n 
is the number of observations. Error statistics were also 
calculated for the models developed by Achim et al. (2006) 
using the current dataset. A fit index (R’) was also 
calculated for each model using the equation given in 
Parresol (1999). Two sets of fit indices were calculated. In 
the first set, the predicted values were estimated from the 
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fixed-effects terms of each model, and in the second, from 
both the fixed and random effects. In addition, plots of 
model residuals versus fitted and explanatory variables 
were visually inspected for any obvious trends. 


3 Results 


Branch diameter ranged from 5 to 73 mm with a mean of 
19 mm, while branch insertion angle ranged from 20° to 
170° with a mean of 78° (Table 2). The mean numbers of 
nodal and internodal branches were 3.7 and 4.9, respec- 
tively. The relative branch diameter was, on average, 0.52+ 
0.28. Overall, there was a significant amount of within- and 
between-tree variation in the branch properties examined 


(Fig. 1). 
3.1 Model evaluation 
3.1.1 Branch number 


The number of branches in an annual shoot increased with 
increasing annual shoot length as well as with increasing 
DBH (Table 3). The number of internodal branches and the 
total number of branches increased with increasing annual 
shoot height in the stem, while the number of nodal branches 
decreased. After accounting for annual shoot length, height in 
the stem, and DBH, re-spacing still had a significant effect on 
the number of both nodal and internodal branches (p<0.05; 
Fig. 2). For nodal branches, the most significant differences 


Table 2 Summary of the branch properties by re-spacing intensity 
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were between the 1.9 and 3.7 m* treatments (p<0.001) and 
the 1.9 and 2.6 m” re-spacing treatments (p<0.05), while the 
4.6 m* re-spacing treatment was significantly different from 
the 1.9 m* re-spacing treatment at the 10% level (p=0.06). 
For internodal branches, the only significant difference 
observed was between the 1.9 and 5.7 m* re-spacing 
treatments (p=0.007). The fixed effects of the models given 
by Eqs. la and 1b explained 30% and 41% of the original 
variation in the number of nodal and internodal branches, 
respectively, and there were no obvious trends in the model 
residuals when plotted against the fitted values or explana- 
tory variables (not shown). 

For the combined number of nodal and internodal 
branches, Eq. Ic explained 42% of the variation in total 
number of branches per annual shoot, while the only 
significant difference observed was between the 1.9 and 
5.7 m* re-spacing treatments (p=<0.001). The mean error 
and mean absolute error calculated from the fixed effects of 
Eq. lc were —0.001 and 2.2 branches per annual shoot, 
respectively (Table 4). In comparison, the corresponding 
model developed by Achim et al. (2006) had a mean error 
and mean absolute error of 1.46 and 3.01 branches per 
annual shoot, respectively, when applied to this dataset 
(Table 4). There was no significant effect of branch status in 
the models for branch frequency. 


3.1.2 Maximum branch diameter 


Maximum branch diameter was positively related to crown 
ratio and had a curvilinear relationship with distance from 


Attribute Re-spacing treatment (m7) 
1.9 2.6 Oe 4.6 a] 

# Observations 530 499 457 515 527 

% Branches alive 41.8 45.4 36.3 44.6 50.1 

# Annual shoots sampled per tree 2.2 12:6 22 128 11.8 

NBy 3.1 (0.9) 3.6 (1.4) 3.94152) 4.6 (1.5) 4.1 (1.1) 
15] [1, 10] [1, 8] (2,9 [2:7] 

NB, 6.5 (4.0) 5.1 (3.4) 4.4 (3.3) 6.2 (5.3) 5.0 (3.7) 
[0, 16] [O, 15] [0, 12] [0, 20] [0, 19] 

BDmax 29 (9.9) 32 (10.1) 29 (9.7) 42 (13.1) 45 (16.0) 
[7, 48] [5, 49] [1, 46] [10, 62] (1, 73] 

BDregr 0.52 (0.30) 0.58 (0.31) 0.59 (0.31) 0.55 (0.31) 0.59 (0.31) 
[0.07, 1.00] [0.02, 1.00] [0.05, 1.00] [0.10, 1.00] [0.05, 1.00] 

VBA 82 (18) 75 (15) 7S (15) 80 (17) 79 (15) 
[30, 170] [25, 150] [25, 140] [20, 140] [20, 140] 


Mean values are presented followed by standard deviations in parentheses and the range of values in square brackets 


NBwy number of nodal branches per annual shoot, NB; number of internodal branches per annual shoot, BDj4x the maximum branch diameter per 
annual shoot (mm), BDpz, relative branch size (BD/BD\ax), and VBA vertical branch insertion angle (°) 
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Fig. 1 Beanplots (Kampstra 
2008) showing the distributions 
of number of branches per 
annual shoot, branch size, 
branch angle, and relative 
branch size by re-spacing 
treatment. Dotted lines overall 
mean, solid lines group means 


Branch diameter (mm) 


Lg 


Branch angle (°) 


the stem apex (Table 3). After these effects were taken into 
account, re-spacing had a significant influence on maxi- 
mum branch diameter. When compared to the 1.9 m? re- 
spacing treatment, maximum branch diameter was signifi- 
cantly larger in the 4.6 m* (p=<0.001) and 5.6 m* (p= 
0.015) treatments (Fig. 3). The fixed effects of the model 
were able to explain almost 75% of the variation in 
BDyjax. Mean error and mean absolute error in BDyax 
were 0.01 and 5.3 mm, respectively, while the mean 
percentage error was 19.7%. There were no trends in the 
model residuals when plotted against the fitted values or 
explanatory variables from the model (not shown), and 
there was no additional effect of branch status on maximum 
branch diameter. 


3.1.3 Relative branch diameter 


Relative branch diameter was positively related to the 
number of branches within an annual shoot and was higher 
in nodal branches, while it decreased with annual shoot 
height in the stem and had a curvilinear relationship with 
branch rank (Table 3). There was no effect of re-spacing on 
relative branch diameter within an annual shoot, and no 
effect of branch status on relative branch size. Overall, the 
fixed effects of the model explained 77% of the variation in 
BDegy, and the mean error, mean absolute error, and mean 
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percentage error were —0.01, 0.11, and 25.4%, respectively. 
There was no trend in the model residuals when plotted 
against the fitted values or explanatory variables (not 
shown). 


3.1.4 Branch insertion angle 


Vertical branch insertion angle had a curvilinear relation- 
ship with distance from the stem apex, increasing rapidly in 
the upper crown before levelling off and then decreasing 
towards the stem base. VBA also decreased with increasing 
branch diameter, the effect of which was greater with 
increasing distance from the stem apex. Insertion angle in 
nodal branches was slightly lower compared with internodal 
branches, while it was generally higher in live rather than in 
dead branches (Table 3; Fig. 4). There was a significant 
difference in branch insertion angle between re-spacing 
treatments after accounting for these effects (p<0.001). In 
general, the smallest branch angles were in the 2.6 and 
3.7 m°* treatments, although only the latter treatment was 
statistically different from the 1.9 m* re-spacing treatment 
(p=0.02). However, the 1.9 m* treatment was not statisti- 
cally different from either the 4.6 or 5.7 m* treatments (p= 
0.822 and 0.992, respectively). Where statistically signifi- 
cant differences were observed, the mean differences were 
generally less than 5°. There were no observed trends in the 


INGA 2 spxnge 





100 


Table 3 Parameter estimates, standard errors and p-values for the 


equations used in this analysis 


Parameter 


NBy (Eq. 1a) 
a) 

a} 

ag 

a3 

a4 

as 

a6 

NB; (Eq. 1b) 
bo 

b 

bo 

bs 

be 

by 

NBr (Eq. Ic) 
Co 

Cj 

C2 

C3 

C7 

BDmax (Eq. 3) 
do 


BDrex (Eq. 4) 
Co 

C] 

C2 

C3 

C4 

C5 

VBA (Eg. 5) 
80 

81 

22 

83 

84 

85 


Y) Springer 


Estimate 


05338 
0.2895 
=0,0095 
0.0130 
O.19T5 
0.2778 
0.2060 


= 10703 
1.6804 
0.0467 
0.0165 

=0,3919 

—0.7403 


0.4828 
0.9338 
0.0190 
0.0160 
—0.4644 


5.8781 
=1.7729 
16.3413 
37,3611 
—1.4772 
=2 1133 
da) 
=2,5019 

0.3553 

2271 

0.8516 

0.8870 


0.8624 
—0.0148 
—0.2808 

0.0186 
—0.0030 

0.1862 


78.4865 
=0.7097 
=1,1167 
10.6581 

0.0216 
—5.3886 





Standard error 


0.1486 
0.0827 
0.0018 
0.0032 
0.0614 
0.0613 
0.0796 


0.3432 
0.1864 
0.0048 
0.0070 
0.1660 
0.2268 


0.1664 
0.0971 
0.0022 
0.0035 
0.1130 


3.8380 
0.1299 
0.6967 
10.9709 
2.0003 
2.0200 
2.6064 
2.5311 
0.1524 
O:1559 
0.1581 
0.1514 


0.0253 
0.0043 
0.0208 
0.0017 
0.0006 
0.0061 


2.2204 
0.0489 
0.1086 
L013 
0.0020 
0.6774 


p-value 


0.0004 
0.0005 
0.0001 
0.0007 
0.0059 
0.0003 
0.0186 


0.0020 
<0.0001 
<0.0001 

0.0311 

0.0297 

0.0043 


0.0040 
<0.0001 
<0.0001 

0.0002 

0.0007 


0.1268 
<0.0001 
<0.0001 

0.0008 

0.4609 

0.2964 

0.0395 

0.3124 

0.0205 

0.1466 
<0.0001 
<0.0001 


<0.0001 

0.0006 
<0.0001 
<0.0001 
<0.0001 
<0.0001 


<0.0001 
<0.0001 
<0.0001 
<0.0001 
<0.0001 
<0.0001 
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Table 3 (continued) 
Parameter Estimate Standard error p-value 
26 5.4260 1.0381 <0.0001 
g7 —4.9653 1.9436 0.0194 


gs =5.9530 1.9569 0.0067 


Nonsignificant (p>0.05) parameters were omitted from the final 
models, except where re-spacing treatment was included as an 
interaction term, in which case all parameters are presented (..e., 
Eq. 3) 


model residuals when plotted against the fitted values or 
fixed-effects variables of the model (not shown). 

The fixed-effects terms in the model developed in this 
study accounted for just less than 33% of the variation in 
branch insertion angle, and the mean error and mean 
absolute error were —0.1° and 8.8°, respectively (Table 4). 
By comparison, the model developed by Achim et al. 
(2006) had a mean error of —2.5° and a mean absolute of 
error of 10.1° when applied to this dataset (Table 4). 


4 Discussion 


Sitka spruce 1s an important commercial species in 
Northern Europe, Canada and the USA, and knowledge of 
its long-term response to different silvicultural treatments is 
needed in order to design effective regimes. The results of 
the current study demonstrate that, in the models for branch 
frequency and maximum branch size, the inclusion of DBH 
and CR, respectively, partially accounted for the effects of 
re-spacing, but significant differences remained, particularly 
in the most extreme treatments (i.e. 1.9 and 5.6 m7’). 
Similarly, Moore et al. (2009) found that branch index 


HB Internode 
Node 





Number of branches per annual shoot 


19m 26m 3.7m 4.6m 5.7m 
Spacing 


Fig. 2 Predicted numbers of nodal and internodal branches (with 
standard errors) per annual shoot using Eqs. la and 1b by re-spacing 
treatment for the mean annual shoot length, annual shoot height, and 
diameter at breast height of the trees sampled in this study. Different 
letters indicate significant differences between groups based on 
Tukey’s test at a = 0.05 
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Table 4 Comparison of mean error (£), mean absolute error (|E]), and mean percentage error (E%) for the equations for branch angle and total 
number of branches developed in this analysis (fixed effects only) and Achim et al. (2006) by re-spacing treatment 


Re-spacing (m7) Branch angle (°) 


Eq. 5 Achim et al. (2006) 
E JE] E% E |E| 
1.9 0.60 0.63 11.70 0.10 10.74 
2.6 —0.49 7.99 10.70 —5.91 9.93 
a7 = 105 8.30 Loy —5.26 10.03 
4.6 —0.96 9.43 L351 O23 10.35 
ja —().64 8.69 10.99 =i 9.37 
Overall —0.10 8.83 25 2.52 10.09 


(defined as the mean of the six largest branches per tree) 
was similar for the three narrowest re-spacing treatments in 
the same re-spacing trial, but was substantially greater in 
the two widest treatments. They also found that the total 
area of knots on sawn timber from this experimental site 
was similar in the 1.9 and 2.6 m* re-spacing treatments, but 
was significantly higher in the 3.7, 4.6, and 5.6 m° 
treatments. 

The presence of a residual effect of spacing in addition 
to that expressed through tree size variables alone indicates 
that the allometric relationships between crown size and 
stem diameter may change in trees grown at wider 
spacings. This has also been observed in other studies 
(e.g. Garber and Maguire 2005; Hein et al. 2008b). In 
practice, this result reveals a potential limitation of a 
modelling approach that links branch characteristics with 
tree growth, as it shows that such relationships may not be 
totally independent of a prior knowledge of the silvicultural 
history of a given stand. This highlights the need to 
calibrate future wood quality models with data from a wide 


Relative height in stem 
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Fig. 3 Maximum branch diameter profiles predicted using Eq. 3 by 
re-spacing treatment for the mean crown ratio values in each treatment 


Total number of branches (nodal+internodal) 


iq. le Achim et al. (2006) 

E% E JE] E% E JE E% 

13.03 —0.001 212 25.43 1.99 3,23 47.99 
12.37 —0.003 2.05 26.52 33 2.58 38.79 
12355 —0.001 2.13 29.02 0.79 2.56 37.63 
12773 —0.001 2.33 22,54 a2 1 3.78 54.43 
35 0.004 241 Die 221 3.05 45.41 
12.44 —0.001 2.20 26.37 1.86 3.01 44.48 


range of silvicultural treatments and stand conditions, and 
illustrates the potential problems of using static assessments 
to represent dynamic processes such as crown development. 

The results of the present study were also generally 
consistent with those of Achim et al. (2006). In both 
studies, branch diameters showed a peak just above the 
crown base regardless of re-spacing treatment, and there 
was no drastic influence of stand density on branch 
insertion angle. Maximum branch diameters were generally 
10-15% larger in the 4.6 and 5.6 m” re-spacing treatments 
when compared to the other treatments. However, in the 
lower crown, the branches could be >120% larger. This 
could clearly affect the visual appearance of sawn timber as 
well as its mechanical properties (Moore et al. 2009). 
Overall, the responsiveness of Sitka spruce branch diameter 
to stand density manipulation was similar to that observed 
in Douglas-fir. For example, using data from a Douglas-fir 
re-spacing trial, Weiskittel et al. (2007) found that there was 
only a ~10% increase in maximum branch diameter in the 
heavy thinning treatment (approximately 250 residual 
stems ha -) relative to the control (~1,000 stems ha °). 

In many studies investigating branch characteristics in 
conifers, the number of branches has generally been one of 
the most difficult attributes to model (e.g. Makinen and 
Colin 1999a, b; Weiskittel et al. 2007), and this was also 
found to be the case for Sitka spruce. The observation that 
the number of branches in an annual shoot was positively 
related to height growth was consistent with the findings of 
Achim et al. (2006) and Cannell (1974). Likewise, Achim 
et al. (2006) found that the number of branches in an annual 
shoot was dependent on its position in the crown even after 
accounting for annual shoot length; however, a curvilinear 
relationship was observed in Achim et al. (2006), while a 
linear relationship was observed in this study. Achim et al. 
(2006) predicted that, on average, the number of both nodal 
and internodal branches increased with thinning, while 
results from this study showed that there was no consistent 
trend between re-spacing intensity and the number of 
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Fig. 4 Predicted vertical branch 
angle profiles for a live nodal 
branches (Eq. la) and b live 
internodal branches (Eq. 1b) by 
re-spacing treatment for the 
mean branch diameter values in 
each treatment 
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branches. Overall, the branch count model of Achim et al. 
(2006) performed reasonably well with the current dataset 
when compared to the models developed in this study, 
despite being developed using data from a more restricted 
range of spacings (1.e., 1.7 to 2.4 m). 

In comparison to other species, Sitka spruce holds many 
branches (Cannell 1974) and often responds to thinning and 
pruning treatments with the production of numerous 
epicormic branches (Deal et al. 2003; Quine 2004). Results 
from this study indicate that creating a widely spaced stand 
with an early re-spacing treatment reduced the number of 
both nodal and internodal branches for a given set of 
covariates, while highly dense stands will have a greater 
number of internodal branches. Hein et al. (2008a, b) found 
that stand density had no effect on the number of branches 
in widely spaced Douglas-fir, while Makinen and Hein 
(2006) found that Norway spruce [Picea abies (L.) Karst.] 
trees growing at 1,600 stems ha’ had slightly fewer 
branches per annual shoot than trees from stand densities 
of 350 and 700 stems ha '. These results indicate that more 
replication might be necessary in order to better understand 
the influence of stand density on the number of branches 
per annual shoot. There could also be a significant 
interaction between treatment and genetics, as Cannell 
(1974) found differences in the number of branches 
between varying seed origins of Sitka spruce. 

Branch insertion angle is an important attribute from a 
wood quality perspective as it determines the area of 
influence on sawn timber for a particular branch. For 
example, a lower branch angle for a given branch diameter 
would result in a larger knot area on sawn timber. The 
models developed in the current study indicated that 
insertion angle decreased with increasing branch diameter, 
which is in agreement with a previous study on radiata pine 
(Pinus radiata D. Don, Grace et al. 1999). However, the 
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authors of the latter study postulated that if longitudinal 
rather than cross-sectional data were available, a positive 
correlation might be expected, a trend which has been 
observed in both Scots pine (Makinen and Colin 1999a, b) 
and Douglas-fir (Hein et al. 2008a). In the present analysis, 
branch angle was significantly influenced by stand density. 
However, the observed differences between re-spacing 
treatments were small (<5°), which is consistent with 
several other studies. For example, Weiskittel et al. (2007) 
found that Douglas-fir branch angle was not influenced by 
intensive management activities, including thinning, and 
Hein et al. (2008b) similarly found no influence of stand 
density in widely spaced Douglas-fir. Despite the high 
observed variation in branch angle in the present study, 
only around 33% of the total variation could be explained 
by the developed model. Regardless, both the model 
developed in this present analysis and that of Achim et al. 
(2006) applied to the current dataset had a relatively similar 
level of performance. 

The distribution of relative branch diameters along the 
stem gives an indication of the range and variability of 
branch sizes within an annual shoot. Re-spacing had no 
influence on the relative diameter distribution within an 
annual shoot. In this study, the distribution of relative 
branch sizes was relatively consistent within a Sitka spruce 
crown, and the decline in relative branch sizes with branch 
rank was steeper than observed for Douglas-fir (Hein et al. 
2008b), probably due to the greater number of epicormic 
and internodal branches in Sitka spruce crowns. 


5 Conclusion 


Overall, this study highlighted the important influence that 
stand density manipulation through early re-spacing can 
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have on branch characteristics, particularly size and 
frequency, which in turn will have an impact on the 
properties of sawn timber. A large part of this variation in 
wood properties was accounted for when a density- 
dependent tree size descriptor (CR or DBH) was included 
in the analysis, but in the case of maximum branch diameter 
and the number of branches per annual shoot there was still 
an additional effect of re-spacing over and above any 
effects attributable to changes in tree-size. This suggests 
that wood quality assessments can be improved if the full 
silvicultural history of the stand is known. It also 
emphasises the need to pursue the development of wood 
quality models by including more data from a wide range of 
silvicultural treatments. The results from this study indicate 
that early re-spacing of Sitka spruce from 1.9 m* to no 
wider than 2.6 m* will not compromise log or sawn timber 
quality in terms of the size and frequency of knots. 
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